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Abstract 

Contemporary attempts to explain the existence of ultra-high energy cosmic rays using plasma-based 
wakefield acceleration deliberately avoid non-Standard Model particle physics. However, such proposals 
exploit some of the most extreme environments in the Universe and it is conceivable that hypothetical 
particles outside the Standard Model have significant implications for the effectiveness of the acceleration 
process. Axions solve the strong CP problem and provide one of the most important candidates for Cold Dark 
Matter, and their potential significance in the present context should not be overlooked. Our analysis of the 
field equations describing a plasma augmented with axions uncovers a dramatic axion-induced suppression 
of the energy gained by a test particle in the wakefield driven by a particle bunch, or an intense pulse of 
electromagnetic radiation, propagating at ultra-relativistic speeds within the strongest magnetic fields in the 
Universe. 


1 Introduction 

The existence of cosmic rays with energies > 10 20 eV remains an outstanding puzzle in astrophysics [l]j2|. The 
well-established Fermi paradigm for cosmic acceleration, in which charged particles gain energy by scattering 
repeatedly from magnetic inhomogeneities [3] or shock waves SHU, has not yet been shown to be efficient for 
accelerating particles to such exceptional energies. However, it has been argued that strategies whose effective¬ 
ness has been honed in the laboratory could be important in the astrophysical context [8]. Experience gained 
in the laboratory suggests that efficient acceleration mechanisms exploit sub-luminal waves, avoid collisional 
processes and do not involve particle trajectories with large curvature (thus avoiding excessive energy losses due 
to inelastic scattering and synchrotron radiation, respectively). 

A sufficiently short and intense laser pulse, or charged particle bunch, propagating rectilinearly through 
a laboratory plasma will drive an inhomogeneity in the plasma electron density that trails behind the pulse 
or bunch. The inhomogeneity, or wakefield, is a wave in the plasma electron density that propagates at the 
group velocity of the pulse (or velocity of the bunch) and whose longitudinal electric field is several orders 
of magnitude greater than that achievable within radio-frequency cavities. Driven mainly by the spectacular 
success of a number of landmark experiments and simulations mm approximately a decade ago, plasma-based 
wakefield acceleration |l3j is now recognised as a vital concept for the next generation of terrestrial particle 
accelerators mm- 

In tandem, plasma wakefield acceleration has also gained attention in astrophysical contexts [ TO fTT ] , For 
example, it has been argued that longitudinal space charge waves excited by Alfven shocks caused by the 
collision of a pair of neutron stars may have the necessary properties for reliably accelerating protons to ZeV 
energies [8J. Thus, plasma wakefield acceleration may provide an effective solution to the problem of the origin 
of ultra-high energy cosmic rays. 

In fact, neutron stars have long been suspects in the search for the culprits behind ultra-high energy cosmic 
rays |18| . In particular, the fast rotation of a pulsar combined with its strong magnetic field yields a particle 
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accelerator based on unipolar induction, and the implications of this mechanism have been thoroughly explored 
over many years (see, e.g., reference ffl)- Although difficult questions remain over whether this approach can 
correctly explain the ultra-high energy cosmic ray spectrum [T] , some success has been achieved in recent years 
in explaining the acceleration of heavy nuclei using sufficiently young strongly-magnetised neutron stars with 
millisecond rotation periods m- Nevertheless, plasma wakefield acceleration is an important alternative to 
unipolar induction, especially for light nuclei mm, and the focus of the present article will be on plasma 
wakefield acceleration in ultra-strong magnetic fields. 

A strength of the plasma wakefield mechanism, like the unipolar inductor, for cosmic acceleration is that it 
does not rely on ingredients outside of the Standard Model of particle physics. However, non-Standard Model 
effects may be relevant as a consequence of the ultra-strong electromagnetic field strengths expected in the most 
extreme astrophysical environments. In particular, it is conceivable that effects due to hypothetical particles 
with very weak coupling to light and matter may manifest in such environments. One of the most popular Dark 
Matter candidates, the axion, was proposed as an elegant solution to the strong CP problem in QCD 12(Tif2T| 
before its significance in the cosmological context was expounded 1331 - 1251 . Furthermore, in addition to the QCD 
axion, light pseudo-scalar particles are a generic consequence of type IIB string theory [26| : hence, a range of 
experiments have been developed, or are under development, in an attempt to uncover the effects of axions and 
axion-like particles (ALPs) 1371 - 1331 . Although positive detection remains elusive in the laboratory, it is possible 
that ALPs play a significant role in the ultra-strong magnetic field of a neutron star. 

The purpose of this article is to demonstrate that ALPs may have significant consequences for plasma 
wakefield acceleration in ultra-strong magnetic fields. Although ALP fields couple directly to ordinary matter, 
the effects of the ALP-photon coupling are expected to dominate over those due to ALP-matter couplings in 
strongly magnetised plasmas. For convenience, we will briefly describe the main ingredients of our approach 
and summarise the key results before turning to their derivation. 


2 Ingredients in the model and key results 

Our assessment of the influence of ALPs is based on a model of a magnetised plasma that includes two charged 
pressureless perfect fluids. One fluid represents mobile electrons whilst the other fluid describes the charge 
carriers of a neutralising background medium with constant proper number density. The charge-to-mass ratio 
of the charge carriers of the background is assumed to be considerably lower than that of the other charged 
particles in the system and we assume that the motion of the background can be neglected over the timescales 
of interest. 

The explanations in references mm of the origin of ultra-high energy cosmic rays invoke acceleration in 
astrophysical jets where the magnetic field is almost certainly too weak for ALPs to play a significant role. 
However, it has also been argued nu that plasma-based acceleration is of interest within neutron stars, where 
the magnetic field is considerably higher. 

The outer crust of a neutron star provides a background medium comprised of magnetically polarised iron 
atoms. Due to the considerable strength of the magnetic field (~ 10 8 T), the electron ‘gas’ within the crust 
is essentially confined to move along magnetic flux tubes threaded by the field lines [16] - Thus, we model the 
plasma wakefield as a sub-luminal non-linear longitudinal plane wave in the electron fluid density propagating at 
velocity v parallel to the field lines of a homogeneous magnetic field of strength B in the frame of the background 
medium. The electron density wave generates an ALP field whose strength is proportional to gB 1 where g is the 
ALP-photon coupling constant; in turn, stress-energy-momentum conservation requires that the plasma fields 
are influenced by the ALP field. A self-consistent analysis of this type was previously used to uncover a novel 
signature of the ALP-photon coupling in waveguide mode spectra [31] ■ 

For sufficiently large density fluctuations, the electric field of the wave has a sawtooth-like profile and its 
amplitude saturates; this leads to an upper bound on the energy that can be gained by a test particle in the 
wakefield. In particular, we will show that a test particle with mass M and charge Q cannot achieve an energy 
greater than W max + Me 2 where 

H 3 g 2 B 2 / tanher 
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with 


Wg^b = M <? [20 2 (7 3 - 7) + 20(7 2 - l)Vl + 0 2 (7 2 - 1) + 7 - l] 
the increase in energy in the absence of the magnetic field. The dimensionless parameters a, 0 are 
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where m a is the ALP mass, 7 = l/i^/l — v 2 /c 2 is the Lorentz factor of the phase speed v of the wave and 
u} p = \J e qono/(£om e ) is the plasma frequency (qo is the charge on a particle of the background medium, no 
is the proper number density of the background medium, e is the elementary charge and m e is the mass of an 
electron). Equation (0 yields the well-established result 

W^ 0 = m e c 2 ( 4 7 3 — 3 7 — 1) (4) 


when the test particle is an electron (M = m e , Q = —e) foolfdlJI . 

The influence of the ALP field emerges as an overall multiplicative factor in dT]) that does not depend on the 
properties of the test particle. It is clear that VF max < Wgl q where the approximate bound is saturated when 
B = 0, and the ALP field and the background magnetic field together reduce the maximum overall energy gain 
of the test particle. The previous result is reassuring because the energy supplied by the driver of the wakefield 
is shared between the ALP field and the plasma, and the test particle does not directly couple to the ALP field 
in our model. Furthermore, as we will show, QED vacuum polarisation does not contribute to ©> to first order 
in small quantities, even though it affects the wavelength and amplitude of the wakefield. 

It is important to note that our analysis does not incorporate the structure and dynamics of the driver of 
the wakefield. Neither the back-reaction of the plasma and ALP field on the driver, nor the back-reaction of 
the accelerated particle on the plasma, ALP field and driver, are included. Detailed investigation of such effects 
requires intensive numerical calculations that are beyond the scope of the present article, but it is likely that 
they would further reduce the energy gain of the accelerated bunch. For example, in most laboratory-based 
configurations, the coupling between the driver and the plasma in particle-beam driven wakefield acceleration 
imposes severe restrictions on the energy gain El- 

The scale of the ALP-induced effects in 0 can be estimated as follows. The DESY ALPS-I (Any Light 
Particle) experiment has excluded light ALPs with strengths g > 10 - 7 GeV -1 , whilst astrophysical and cos¬ 
mological considerations lead to 10 -6 eV < m a c 2 < 10 -2 eV for the QCD axion m- The plasma frequency 
of the magnetically polarised iron lattice providing the neutralising background within the outer crust of a 
neutron star is ~ 10 18 Hz when B ~ 10 8 T m ■ Hence, convenient estimates of the dimensionless parameters 
that characterise the ALP-induced suppression of a plasma wakefield accelerator driven by a ~ 1 TeV bunch of 
electrons (7 ~ 2 x 10 6 ) propagating along the ~ 10 8 T magnetic field within the outer crust of a neutron star 
are as follows: 


h 3 g 2 B 2 
Mo c m 2 


cr = 9.7 
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Note that a ~ 3% perturbation to IFjf a Q follows from ©, 0, @ when the representative parameters are 
used. The relative size of the perturbation is highly sensitive to the ALP mass; in particular, it reduces to 
~ 0.04% for m a c 2 = 10 -4 eV and increases to ~ 90% for m a c 2 = 10 - 6 eV. As we will see in the remainder 
of this article, a perturbative approach is used to obtain 0 and, although the use of perturbation theory is 
suspect when m a c 2 < 10 -6 eV, it is reasonable to conclude that the ALP field has a substantial effect on the 
maximum energy gain if m a c 2 < 10 - 5 eV. 

More stringent upper bounds on the ALP-photon coupling constant 3 , such as those obtained from solar axion 
searches by CAST [38] (g < 10 -10 GeV -1 ), can be compensated by the magnetic fields found in magnetars (B ~ 
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10 11 T). The plasma frequency of the outer crust of a neutron star satisfies lo p ~ 2n x 10 18 (B/10 8 T )‘ 5 rads -1 , 
where 5 ~ 3/5 — 3/4 [T 6 ], and inspection of © shows that a wakefield driven by a ~ 30TeV bunch of electrons 
in a ~ 10 11 T field leads to values of a that are similar to those found when B ~ 10 8 T. Note that only 
relatively modest drive-bunch energies (in the context of cosmic rays) are required to obtain substantial axionic 
suppression of the plasma wakefield acceleration mechanism. 

3 Derivation of the results 

Henceforth, units are used in which the speed of light c, the permittivity £o of free space and the reduced 
Planck constant h are unity. Furthermore, for linguistic convenience, we will refer to the charge carriers of the 
neutralising background medium as ions regardless of whether they are polarised ion cores in the crust of a 
neutron star or ions in a magnetised plasma. 

3.1 Field equations 

The variables describing a cold ALP-plasma are the ALP 0-form a , the electromagnetic 2-form F , the 4-velocity 
field V e of the plasma electrons, the proper number density n e of the plasma electrons, the 4-velocity field Vo 
of the ions and the proper number density no of the ions. The effect of the electromagnetic field on the ions is 
negligible over the length and time scales of interest; we choose no to be constant and choose Vo = d/dt where 
the spacetime metric r/ and volume 4-form *1 are 

rj = —dt (g> dt + dx ® dx + dy g) dy + dz ® dz, (7) 

* 1 — dt /\dx Ady Adz ( 8 ) 

with A the exterior product and d the exterior derivative on differential forms. Exterior differential calculus is 
used extensively in this section because it is an efficient tool for formulating the ALP-plasma field equations 
and reducing them to non-linear ODEs; a detailed account of the techniques and conventions used here may be 
found in reference {551 . 

The Hodge map * induced from the volume 4-form *1 satisfies the identity *(/? A W) = iw * /3 for all 
differential forms /3 and vectors W. The action of * is extended to non-decomposable forms by linearity, and 
the linear_operator lw is the interior derivative with respect to W. The 1 -form W is the metric dual of W and 
satisfies W(V) = y(W, V) for all vectors V; likewise, the vector field f3 is the metric dual of the 1 -form /? and 
satisfies V(fl) = /3(V) for all vectors V. 

The electric 4-current densities of the plasma electrons and the ion background are q e n e V e , qonoVo , respec¬ 
tively, where q e = — e is the charge on an electron, qo is the charge on a background ion and r](V e ,V e ) — 
rj(V o, Vo) = —1 . The behaviour of the electromagnetic field is determined by the Gauss-Faraday and Gauss- 
Ampere laws: 


dF = 0, (9) 

d-k G — -q e n e * V e — qono * Vo (10) 

and the ALP field a satisfies 

d-k da — mr a a * 1 = — d a X * 1 (11) 

with m a the ALP mass. The electromagnetic excitation 2-form G is specified by 

G = 2(d x XF-dyX * F) (12) 

where A is a 0-form-valued function of the ALP field a and the electromagnetic invariants 

X=k(FAkF), Y = k(F A F). (13) 
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The O-form A(A1, Y, a) is the sum of all purely electromagnetic contributions to the Lagrangian (including 
effective self-couplings due to QED vacuum polarisation) and terms that encode the interaction between the 
electromagnetic and ALP fields. 

The behaviour of the electron fluid is determined by appealing to total stress-energy-momentum conservation. 
The divergence of the total stress-energy-momentum tensor of the electron fluid, electromagnetic held and ALP 
held must balance the forces on the ions. For present purposes, it is useful to cast this statement in the language 
of differential forms; the behaviour of the plasma electrons is determined by the following held equation: 

dr K = q 0 n 0 L K F A *V 0 (14) 

where the total stress 3-form tk of the electromagnetic held, ALP held and plasma electron fluid is 

tk = ikF A *G + A * K + i (i^eda A *da + da A lk * da — m 2 a a 2 * K ) + m e n e rj(V e , K)*V e (15) 

with K a Killing vector and m e the mass of an electron. The stress 3-form tk is related to the total stress- 
energy-momentum tensor T of the electromagnetic held, ALP held and plasma electron huid via the identity 
F{K, W) = L\v * t k, for all vector helds W, and (fTTI) is the component of a local balance law associated with 
K . In particular, if if is a generator of spatial translations then ed is the if-component of a held equation 
describing the local balance of linear momentum whilst a generator of time translations leads to local energy 
balance. 

Before reducing the above held equations to a system of non-linear ODEs, it is worth commenting on the 
approach adopted here in comparison to our earlier analysis of plasma wakehelds in non-linear electrodynam¬ 
ics gD]. It can be shown that the system ©, (flOl) . (fill) , (fill) is equivalent to ©, (ITU) . (fTTI) with the Lorentz 
equation of motion for the plasma electron huid: 

Ly e dV e = —iv e F. (16) 

m e 

The system ©, (USD, ED is the starting point for the analysis in reference jTUj; however, the analysis in 
reference [TO] includes the derivation of a hrst integral that can be shown to follow immediately from the stress- 
energy-momentum balance law (ED EH- For present purposes, it is more efficient to adopt the strategy that we 
recently developed in reference mi and begin with the stress-energy-momentum balance law ED rather than 

ED- 


3.2 ODE system describing non-linear longitudinal ALP-plasma waves 

Solutions to ©, ED. ED. ED are sought that describe longitudinal plane waves propagating along a constant 
magnetic held. The electromagnetic held F and electron huid 4-velocity V e have the form 

F = E((/) dt Adz — B dx A dy , (17) 

K = MC) 9 1 — \J y(C) 2 — 7 2 0 2 (18) 


where ( = z — vt with v the phase speed of the wave, 0 < v < 1 , 7 = l/\/l — v 2 , 6 1 = —dt + vdz and 
9 2 = d£ = dz — vdt. Likewise, the electron proper number density n e and the ALP held a are assumed to 
depend on £ only. 

Clearly, the Gauss-Faraday law m is trivially satished by ED- Furthermore, the point-wise behaviour of 
the components of * G with respect to the co-frame {dt, dx, dy, dz} depends on (( only and it immediately follows 
that dC, A d-kG = 0. Thus, the exterior product of dC, and the Gauss-Ampere law (flOl) yields 


n e 


yo 2 
-nouq 

9e 


1 


(19) 


The set of Killing vectors {d/dx, d/dy, 9 1 , 0 2 } is a basis for vector helds on hat spacetime. Setting K to d/dx 
or d/dy in (fl5l) yields tk — 0 where ~ denotes equality modulo exact forms. On the other hand, inspection of 
ED reveals that the right-hand side of (fTTI) vanishes when K e {d/dx, d/dy} and so (fTTp is trivially satished. 
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Setting K = 6 1 in (fl5l) reveals 
tk — —m e —novy 6 1 A dx A dy 

Qe 

where m has been used to eliminate n e , and 

„ TOp 1 , 

E = —~V 

Qe 7 

follows from ED, (ED, ED- a prime indicates differentiation with respect to 

The remaining outcome of ED is revealed by choosing K = 9 2 . It can be shown that 


( 20 ) 


( 21 ) 


t k ^ 


2(—<9 y A E 2 - dy A EB) + A + 


1 a' 


- m n a - 


qo 

<7e 


m e n 0 v^/y 2 - y 2 


8 1 A dx A dy 


( 22 ) 


where (fl9l) has been used to eliminate n e . Thus, combining (E2l) with the result of using (12T1) to eliminate E in 
the right-hand side of ED reveals 


0 = 


2(<9 a -A E 2 + d Y A EB) — A - 


1 fa' 


*1 2 
— m„a 


+ —m e n 0 (vy/y, 2 - y 2 - y) 

Qe 


(23) 


Only the ALP field equation dill) remains and it can be shown that 


= -8 n X. 
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(24) 


It is worth noting that the above procedure bypasses the second-order ODE for y that emanates from the 
Gauss-Ampere law ED- However, the latter ODE also follows as a consequence of ED, ED. ED and no new 
information is uncovered. 

In summary, the ALP-plasma field system reduces to ED. ED with E given by ED- The electromag¬ 
netic invariants are X = E 2 — B 2 , Y = 2 EB. For definiteness, from now on we will only consider theories 
with a minimal coupling between the ALP and the electromagnetic field rather than those with more general 
couplings (3TJI121. We will focus on theories of the form 

(25) 


A(X, Y, a ) = A E m(AT, Y) + \jaY 


where Aem(A1, Y) depends only on the electromagnetic field invariants and g is a coupling constant. The term 
gaY/2 coupling the ALP field and electromagnetic field does not appear in the stress-energy-momentum tensor 
and its contributions to the first three terms in (1231) cancel; substituting A in (1^31) . (EH) leads to 


0 = 


2(<9aAem E 2 + cVAem EB) - Aem — 


1 ( a 1 


12 


^2 2 
- m„a 


+ —m e n 0 (v^/y 2 - 7 2 - y) 


(26) 


and 


^ - m 2 a a = -gEB. (27) 
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The 0-form Aem (AT, Y) reduces to the classical vacuum Maxwell Lagrangian X/2 in the weak-field limit and, for 
stronger fields, captures the effects of quantum vacuum polarisation. As we will soon see, the detailed structure 
of Aem (AT, Y) is unimportant for our purposes. 
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3.3 Periodic solutions 


Inspection of (071) shows that, in general, a behaves exponentially when g = 0 and such solutions describe ALP 
fields generated by a source other than the plasma electrons and ions. Although such effects could be attributed 
to the driver of the plasma wakefield, we will only consider the ALP fields self-consistently generated by the 
plasma electrons and ions. 

The structure of (l26l) suggests the existence of solutions for g that are periodic in ( when a = 0. Furthermore, 
the structure of (071) ensures that an ALP held generated by a periodic electric held will be periodic; hence, in 
the absence of a background ALP held, a Fourier series representation for g, a will be sought: 


OO 

m( 0 = Vn exp 

n =—oo 

OO 

a(C) = ^2 a n exp 

n =—oo 



(28) 

(29) 


where the period l of the solution must be determined as part of the analysis. Thus ED, ED, E3, E3 lead 
to the relationship 




gBl m e 
47 r 2 n 2 + m 2 7 2 / 2 q e 


2ning n 


(30) 


between the Fourier components of a and g, and we can express (1261) as the non-linear and non-local equation 


0 = 


2(<9vAem E 2 + cVAem EB) 


Aem + <p[f .<] + —m e n 0 (v^/g 2 - q 2 


<7e 



(31) 


for g where the functional ip is 

( 32 ) 

with a specihed by (1711) . (PHI) and 




exp(—27 rin(/l) g(Q d(. 


(33) 


3.4 Maximum energy gain 

Following the approach used in reference EB. the change in energy A£k of a test particle between the endpoints 
I, II of a segment C of its world-line in an inertial frame of reference adapted to a timelike unit Killing vector 
K is 

A£ k = -Mr / (C,K)\ 1 1 1 

= Q [ lkF (34) 

Jc 

where M, Q are the mass and charge of the test particle, respectively, and C is the 4-velocity of the particle. 
Equation (1371) is a covariant expression of the relationship between the change in energy of the particle and the 
work done on the particle by the Lorentz force between I and II. Of particular interest here is the case when 
Ci, Cii are located at adjacent turning points of g, i.e. adjacent nodes of the electric held. 

The square root in (1311) ensures that g > 7 and leads to an upper bound on the amplitude of the solution 
to (15T1) if Aem = X/2 and g = 0 (i.e. the effects of quantum vacuum polarisation and the ALP are neglected). 
However, the effects of quantum vacuum polarisation and the ALP are expected to be small and, from a 
perturbative perspective, an upper bound on g is also expected to exist when such effects are included. 
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The electric field lkF in a frame adapted to K = 7 (dt + vd z ) (i.e. a frame in which the wave is static) is 
E"/d£ using m and dQ = dz — vdt. Thus, the maximum gain in energy of the test particle is 


A £k 


Qme /^d/x 

7 <Ze J Cl 


Qrn e 

- {m - H 1 ) 

7 Qe 


(35) 


where m has been used. Although (1351) can be expressed in terms of the potential difference between two 
adjacent nodes of the electric field, it is convenient for present purposes to retain A Ek in the form given above. 

The parameters /q, /in are adjacent turning points in /i and further analysis is required to determine them. 
Integrating m immediately yields 


|</)[/x] + —m e n 0 (v^n 2 - y 2 

l Qe 



(36) 


because E oc /x' and so -E(Ci) = -E(Cn) = 0. In order to ensure that the test particle is accelerated, rather than 
decelerated, we require the particle to be in the appropriate phase of the wave. If Q/q e > 0 then we choose /q 
to be a minimum of /i and /in to be the subsequent maximum of /i. However, if Q/q e < 0 then we choose /q to 
be a maximum of /i and /in to be the subsequent minimum of /i. 

Note that non-linear corrections to Aem = A /2 only contribute to (l36l) through However, as noted 

above, the corrections to classical Maxwell theory due to quantum vacuum polarisation and the coupling to the 
ALP are expected to be amenable to perturbative analysis. Since </>[/x] is a small perturbation to the remaining 
terms in (1361) , we expect the value of /i for a classical cold plasma to be adequate for calculating the first term in 
(1361) . Furthermore, inspection of (1361) reveals that the absolute minimum value of /x is 7 when </>[/x] is neglected 
and this conclusion will hold when </>[/x] is not neglected. Thus, we express (1361) as 


pMI? 1 + e—m e n 0 (vy/n 2 - 7 2 - /i» + 7) « 0 

^ Qe 

with v the solution to (1261) for /x that arises when Aem = A/ 2 , a = 0. The parameters e, //,* satisfy 


(37) 


e 


+1 if Q/q e > 0 
-1 if Q/q e < 0 


and 


= 



if Q/q e > 0 
if Q/q e < 0 


with /q = 7 if Q/q e > 0 or /xq = 7 if Q/q e < 0 . Setting Aem = A/ 2 , a = 0 in (l26l) reveals 



-7 4 w^(uvV 2 -7 2 



= 0 


(38) 


(39) 


(40) 


where lo p = y/—qoq e no/m e is the plasma (angular) frequency, and it follows that v satisfies 

- 7 4 w 2 (u-\/ v 2 - y 2 - v + 7 ) = 0 (41) 

since, by construction, v' = 0 when v = 7 . 

Note that in calculating the maximum gain in energy, all electromagnetic self-couplings (such as those that 
arise from QED) are distilled away without further approximation. The explicit modifications to the period 
and profile of the wave due to vacuum polarisation conspire to produce no change in the maximum possible 










energy gain [4Tj. The effects of quantum vacuum polarisation implicit in the remnants of the ALP field in (l36ll 
are lost when is approximated by ip\v]. 

Solving (1371) for /i* and casting the result in an amenable form is a little involved and, to avoid distraction, 
we have presented the details in the Appendix. We find 


g* 


s=o 


= 7 3 (1 + v 2 ) 


in the absence of the ALP-photon coupling and 
2 g 2 B 2 


M* « 7 d (l +v 2 ) - 


-7 


2 , 

^ IT s\ 

1-tanh 


7TS 

V 2 /. 


to first order in g 2 , where the parameter s is 

rria'yl 
S = 2 tt 


(42) 


(43) 


(44) 


and only the leading order 7 -dependence of the 0(g 2 ) term in (l43l) has been retained. Hence, using (1521) and 
(1551) with gi = 7 if Q/q e > 0 or gu = 7 if Q/q e < 0 , the maximum gain in energy in the wave frame when g = 0 
is 


A £k 


3=0 


2m e 7 2 n 2 


Q 

9e 


and 


A £k 


2m e 7 2 


Q 

Qe 



g 2 B 2 

m 2 a 



(45) 


(46) 


follows from a perturbative analysis of (143[) . (1551) in g 2 . 

The maximum energy gain in the inertial frame in which the ions are at rest is straightforward to obtain 
by adapting the approach in reference m to accommodate a test particle with arbitrary mass M and charge 
Q. Unlike (|46[) . the energy gain in the ion frame depends on the initial velocity of the test particle. Although 
it is intuitively obvious that the test particle should begin at rest in the wave frame to gain maximum energy 
within the wave, we will prove this result for completeness. 

Expressing the 4-velocity C of the test particle as C = T(K + UL) where (0,0, XJ) is the 3-velocity of the 
test particle in the wave frame, K = ^{d t + vd~), L = y(9 z + vd t ) and T = l/Vl — U 2 gives 

A£ Bt =-Mr t (C,d t )\? 

= M 7 [r„(l + U n v) - ri(l + Uiv)] (47) 


where C/i, U\\ are the values of U and T i, Tn are the values of T at the spacetime events I, II respectively. The 
expression 


A£ k = -Mt,(C,K)\Y 
= M(r„ - Ti) 


(48) 


fixes Tu (and therefore also fixes U\\) in terms of Ti and the known quantity A £k- Since A £k is independent 
of Ti it follows that dTn/dri = 1 using (1551) and 


dA£ dt 

dl'r 




(49) 


follows from (1571) : therefore dA£Q t /dTi < 0 since U\\ > U\ (the test particle is assumed to be in the phase of the 
plasma wave in which it is accelerated). Hence, the maximum value A£]^ ax of A£g t is obtained at the minimum 
value of Pi and setting Ti = 1 in (1571) yields 


A£%™ = M 7 [r n (l + U llV ) - 1] 


(50) 
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( 51 ) 


where Tn is specified in terms of A £k as 

Fn = 1 + Ti AE k - 

The maximum total energy that the test particle can achieve in the ion frame is TU max + M where 

^rmax = A £max + M ^ ^ (52) 

Using (145II , (l50ll , (IFTl) to eliminate A£JJJ ax from (1521) in favour of 7 and the ALP-plasma parameters leads to 
lT max | s=0 = M[ 20 2 ( 7 3 - 7 ) + 20( 7 2 - 1) y/l + 0 2 (7 2 - 1) + 7 - l] (53) 

where 


0 = 


m e 

M 


and a perturbative analysis of (15(1 in g 2 yields 

max ^ j^max 


1 J 

L 9 2 b 2 \ 

0 

II 

©5 

l rnl [ 


, 2 1 , fns\ 

1-tanh — ) 

7 r s 12/ 


(54) 


(55) 


for 7 1. The equality in (1551) is approximate because it is only valid to first order in g 2 and the dominant 

behaviour of the 0(g 2 ) term in large 7 has been exploited to write the final result in a factored form. Equation 
(El with A/ — /y/.f.. (J — (](: — —e leads immediately to the classic result 


W n 


ls=o 


= m e (47 3 - 37 -I) 


for the maximum energy gain of a test electron in a plasma wakefield accelerator 
in the Appendix, it follows that ns/2 m a where 


V2 


m a 7' 


.3/2 


and hence 

VU max 


UJri 




(56) 

]. Using El and (1641) 

(57) 


I g=0 


1- 


g 2 B 2 


1 - 


tanh er 


(58) 


The starting point in Section d is obtained by restoring appropriate powers of c, ^o, h in El i El i El - 


4 Conclusion 

Axions are hypothetical particles that solve the strong CP problem and are one of the most promising candidates 
for Cold Dark Matter; it is conceivable that they pervade the cosmos. In addition to QCD axions, type IIB 
string theory predicts a plethora of axion-like particles (ALPs) and we have unearthed some of the implications 
of ALPs for plasma-based wakefield acceleration in ultra-strong magnetic fields. Typical values of the ALP 
mass and ALP-photon coupling strength suggest that minimally-coupled ALPs generated within the plasma 
wave could dramatically suppress the effectiveness of the acceleration process in ultra-strong field environments. 
While such suppression is unlikely to play a significant role in scenarios that exploit moderate magnetic fields, 
such as those in references mm, the theoretical significance of axions, and ALPs in general, suggests that 
such particles should be taken into account in any attempt to explain ultra-high energy cosmic rays using 
plasma-based wakefield acceleration in the strongest magnetic fields in the Universe. 

Our analysis does not include the evolution of the the driver at the head of the wakefield, ALPs directly gen¬ 
erated by the electromagnetic field of the driver, or the ALPs generated by external sources; such developments 
could form the basis of a future study. 
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Appendix: The value of /i* 


The maximum amplitude electrostatic wave in the absence of quantum vacuum polarisation and the ALP field 
is governed by 

^ = \f 2 ') 2 u) p {vy/v 2 - q 2 - v + 7) 1 / 2 ( 59 ) 

ac, 

where the electric field is given by CD with fi replaced by v. It is straightforward to obtain an approximate 
analytical solution to (l 59 l) and determine a good approximation to the Fourier coefficients a n of the ALP field 
when 7>1. Inspection of CD reveals that the maximum and minimum values of v are 7 3 (l+u 2 ), 7 respectively 
and, in the following, we will choose the constant of integration in ( 15 U 1 ) such that v = 7 at £ = 0. 

The phase (] of the wave can be written as 


CM 


_jl_ r _ i _ dx 

V / 27 2 Wp Jj (vy/x 2 - 7 2 - X + 7) 1 / 2 



( 60 ) 


where v = z^/7 3 , % = x/7 3 and v = \Jl — 7 2 has been used to clearly show the 7 dependence of the integrand. 
However, 



0 { 7 ” 4 ) 


( 61 ) 


12 












and hence 



is the dominant behaviour of (l60ll when 7 > 1. The period l of the oscillation is 


i=2cr j3{1+v2) 


and thus 


l 


2 cr =2 

s \p=0 


V 27 

UJp 


follows from (1671) . Solving (1()71) for v and writing the result in terms of l using (1641) yields 


KC) 



for 0 < ( < l 


when 7 1 . 

Thus, the coefficients v n of the Fourier series 


( 62 ) 


(63) 


(64) 


(65) 


OO 

KC) = Vn exp 

n=—oo 



follow immediately: 



2nin^/l) i/(£) d(' 

for n = 0 
for n/ 0 


( 66 ) 


(67) 

( 68 ) 


when 7 1 . 

Although (1651) is a poor approximation to v very close to the minima of 17 numerical investigation reveals 
that the difference between the approximate value of (1671) and the exact value is less than 1% for |n| £ { 0 , 1 , 2 ,3} 
with 7 = 10. The percentage error slowly increases with \n\, but the Fourier coefficient v n rapidly decays with n 
and the absolute error is negligible. Moreover, the accuracy of the approximation increases with 7 ; laser-driven 
plasma wakehelds have a Lorentz factor 7 in the range 10 — 100 and the Lorentz factor of laboratory-based 
electron-driven plasma wakehelds are even higher (7 ~ 10 5 ). For notational convenience we will henceforth 
treat (l68l) as an equality rather than approximate equality. 

We now turn to the value of the remnant ip[v\ of the ALP Held in (1371) . Inspection of (1^81) and the expression 


Q-n — 


gBl 


m e 


27 t invri 


47 r 2 n 2 + m^ 7 2 / 2 q e 
obtained from (1301) reveals that = —a n and so, using (1751) . 


(69) 


OO 

a(0) = ^2 a n = 0, 

n =—00 
00 

a(l/2)= J2 M-l) n = 0 . 

n =—00 


(70) 

(71) 
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It is also straightforward to show that 


... v—-v 2nin 

a (0) = 2^ a n —— = K e + n 0 , 

TI— — 00 

a'(l/ 2)= £ o n -j^(-l) n = «e-« ( 


(72) 

(73) 


where K e , result from splitting the Fourier sums into even and odd indices, respectively: 


E 

n even, n ^ 0 

m e 4y 3 \ - 1 


1 


l tie up \ ' 
^o=gB-- T ^ 


g e 7 r^ z —' n 2 + s 2 

n odd 


The details of ED, ED follow from (ED , (ED with the parameter s given as 
m a jl 


s = 


2-7T 


Furthermore, the summations in (1741) . (1751) can be expressed in closed form: 
1 7r coth(7rs/2)s — 2 


E 

n even, n ^ 0 

E ~ 2 ~ 

Z—/ 4. 


U 2 + S 2 


2s 2 


odd 


n 2 + s z 


7r tanh(7rs/2) 
2s ‘ 


Evaluating (1571) at g, = v and using (1701) . (1711) . (1771) . (1731) leads to 
l Cn j 2 «: e K 0 / 7 2 if Q/q e > 0 


if[i 


Cl l- 2 /t e K 0 / 7 2 if < 5 /g e < 0 


where 


and 


Ci = 


Cn = 


0 if Q/q e > 0 

1/2 if Q/q e < 0 


1/2 if Q/q e > 0 
l if Q/q e < 0 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 


have been used, and the periodicity of the fields has been exploited. 

Finally, noting that n e K 0 = 0(g 2 ), equation (1731) follows from a perturbative analysis of (1371) . Equations 

ED, ED yield 




7 3 (1 + v 2 ) - 


9 9 


(82) 


to lowest order in the ALP-photon coupling constant g , and substituting (1771) . (1751) . (1771) . (1751) in (1571) leads to 
(1751) . 

Although only the leading order dependence on 7 has been retained in the 0(g 2 ) term in (1571) . we have 
retained the exact 7 -dependence of the remaining terms for accuracy. 
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